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We clarify general properties of the energy gap regarding its functional dependence on the energy- 
momentum dictated by the invariance under a space inversion or a time reversal. Then we derive 
perturbatively the equation of the imaginary part of the gap function for dense QCD in weak 
coupling and generalize our results from 2SC case to CFL case. We confirm that the imaginary part 
is down by g relative to the real part in weak coupling. The numerical results show that, up to the 
leading order, the imaginary part is no larger than one MeV at extremely large densities and can 
be as large as several MeV for the densities of physical interest. 
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(D ■ I. INTRODUCTION 

\ At sufficiently high density a quark matter is expected to become a color superconductor, which has stirred a lot of 
, interests The mechanism of color superconductivity is essentially the quark analogue of BCS scenario^, which 

implies that if there is an attractive interaction in a cold Fermi sea, the system is unstable against the formation of 
^1 : a diquark(Cooper pair) condensate 0, 0, In QCD case at asymptotically high density, the dominant interaction 
ly-j between two quarks is due to the one-gluon exchange, which is attractive in the color antitriplet channel. A first 
^ principle calculation can be carried out perturbativelyat this den sity and the asymptotic formula for the energy gap 
■ and the transition temperature have been derived 0,11, 0, 0, 0, 0, 0, Q, 0| . At moderately high density that is 
\ accessible in a compact star, the instanton mediated interaction may dominate, which is again attrarctive in the color 
i antisymmetric channel. Because of its nonperturbative origin, one has to resort to various phcnomcnological models, 
' such as NJL model [T^ ITtI ITsL IT^ to explore the color superconductivity in this case. 

. Color superconductivity has a very rich phase structure because of its color and flavor as well as spin degrees of 
• freedomQ. It is well-established that at extremely high density, where the chemical potential is much larger than the 
I strange quark mass, the ground state is characterized by the so-called color-flavor locked(CFL)condensate0. But 
Q ■ actually the chemical potential is unlikely to be much larger than 500Mev in the cores of a compact star. Therefore, the 
^ I heavier strange quark may not be able to participate the pairing with up and down quarks. The color superconducting 
' phase with only two light quarks is normally called the 2SC phase. However, because of the /3-equilibrium and the 
' charge neutrality conditions of the quark system, nontrivial relations will be imposed between the chemical potentials 
• ^ of different quark flavors '2 1]. In turn, such relations could substantially influence the pairing dynamics. It was shown 
that gapless 2SC 22, 23, 24jand gapless CFL|i2^]as well as gapless LOFF|26i|could appear in quark matter regarding 
^ the influence of /3-equilibrium and the charge neutrality. 
■ - - I One of the main consequences of color superconductivity in dense quark matter is the opening of an energy gap, 
which impacts on the transport and thermodynamic properties of quark matter. The general pairing potential is 
retarded and contains damping terms (Landau damping in case of one-gluon exchange of QCD), which corresponds to 
branch cuts along the real axis of the complex energy plane. Consequently, the gap function depends on the energy- 
momentum and acquires a nontrivial imaginary part along the axis of real energy. A gap function with a nonzero 
imaginary part is actually well known from strongly coupled electronic superconductors, as studied in Eliashberg 
theory[23. But the parallel case has not been analyzed in detail until recently by Rueter [23 • Our main work in this 
paper is to derive the equation of the imaginary part of the gap function in weak coupling. Our results confirm that 
the imaginary part of the gap function is down by g relative to the real part in weak coupling. At extremely large 
densities, the imaginary part of the gap function is smaller than one MeV. For chemical potential that are of physical 
interest, /i <1 GcV, the imaginary part can be as large as several MeV. 
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The current work is organized as follows. In section 2, we shall demonstrate some general properties of the energy 
gap as a function of the energy and momentum. Then we review the Eliashberg equations for the energy gap at real 
energy in Section 3. Its solution will be presented in Section 4 and the concluding remarks will then be given in 
Section 5. Moreover, technical details on the derivation of the real part of the gap will be sketched in Appendix A. 
Our units are ?i = c = fc^ = 1 and 4-vectors are denoted by capital letters, K = — (fc, uj) in our formulas. 

II. GENERAL PROPERTIES OF THE ENERGY GAP AS A FUNCTION OF ENERGY AND 

MOMENTUM 

Before embarking on the calculation of the imaginary part of the energy gap, we shall clarify some general properties 
of the energy gap regarding its functional dependence on the energy and momentum dictated by the invariance under 
a space inversion, 7^, or a time reversal, T. These properties are model independent and exact to all orders of the 
coupling strength. While, some of them may be known in the case of electronic superconductivity, they remain murky 
for color superconductivity. Since a large body of formula employed in this section stems from that of Ref . [2fll| . we shall 
follow the convention there for the gamma matrices (i.e. all gamma matrices are hermitian with 71, 73 antisymmetric 
and 72, 74 symmetric). 

Starting with the Schroedinger representation, the Nambu-Gorkov(NG) quark field reads 

and is self-conjugate, i.e. 

*(r,0) = C^'^(r,0), (2) 
where "0 is a Dirac spinor, ipc = 172^'^ with the tilde standing for transpose and 

C = ^^l ^^)=C = -CK (3) 

Color and flavor indexes are suppressed. We choose the phase of ijj such that the transformation rules under V and 
T are given by 

V^{r, 0)r^ = Up^{-f, 0), (4) 

V'i\r,0)V^ = «'^(-f,0)t/^, (5) 

T^{f,0)T-^ = t/T^'(r,0), (6) 

and 

T^\f,Q)T-^ = *^(f,0)C/j,, (7) 

where Up = diag.(74, — 74) — Up, Ut — diag(i777i 73, 177*7173) — —IJt with 77 a phase factor {\ri\ = 1) independent of 
the coordinates and internal indexes. 

The self-conjugation Q and the transformation rules Q-Q in the Heisenberg and Matsubara representations 
follow from the real(imaginary) time development formula of an operator O, i.e. 

0{t) = e'"*OiO)e-'"' (8) 

and 

0(r) = e"^O{0)e~"^ (9) 

with H the Hamiltonian. Note that the similarity transformation (jSJ is not unitary and thus does not preserve the 
relation of hermitian conjugation. In particular, on writing 

6{T) = e"^O\0)e-"^ (10) 

we have 0(r) — 0^{~t) and vice versa. The over bar in this section does not mean a right multiplication by 74. 



A. Real time formulation 

The retarded NG quark propagator is defined by 

S^p{r, t)^i< {vl/„(f, i), -^{{Q, 0)} > {Ti)),pe{t) (11) 

where r4 = diag(74,74) and < O > denotes the thermal average of the operator O, i.e. 

Tre-'3^0 , , 

< O >^ ^^-^ (12) 

with (3 = 1/T the inverse temperature. The Fourier transformation of (|ll(l reads 

S{p,uj)= f dt ( d^re"^'"'P-^S{r,t) (13) 



with (p, Lo) the momentum and the energy. On writing 

The off-diagonal matrix elements carry the long range order and the gap functions for 2SC are extracted according 

= -'7*'/'(P,^)75A2T-2 (15) 

and 

*'(P, t^) = ?70'(P, w)75A2r2 (16) 

with A2 the second Gell-Mann matrix acting on colors and T2 the second Pauli matrix acting on flavors. 
The self-conjugate relation (j^J in the Heisenbcrg representation reads 

^{r,t) = c¥{r,t). (17) 

Taking the complex conjugate of Hll|l and applying the relation H17() we find 

Slp{r,t) = -*CAp(Ct)^„ < |{M'^(f,t),vI/t(0,0)}| > (r4)A;3 (18) 

Suppressing the indexes and using the relations C — C and r4Cr4 — — C, we obtain that 

S*{r,t) = C^S{r,t)C (19) 

and its Fourier transform 

S* i~p, -oj) = C^S{p, uj)C. (20) 
The same relation should also apply to the inverse propagator, i.e. 

S*-\-p,-oj) = C^S-'^{p,uj)C, (21) 

which implies that 

0'(p,c^) = r(-P,-c^). (22) 

Therefore, ^{p,uj) — -~ri*(f){p,uj)j5X2T2 and $'(p, w) — 7]4>* {—p, —uj)j5X2T2 We are left with only one complex gap 
function, <j){p,uj), to consider. 

The V and T transformation rules of the Heisenbcrg operators ^(r, t) and '^^{r,t) read 

V^{r, t)V^ = Up^{~r, i), (23) 
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V¥{f,t)r'' = ¥{-'r,t)Ul, (24) 

r*(f,i)r-i = [/T«'(f,-i), (25) 

and 

r*^(f,i)r-i = ■^\r,-t)ul. (26) 
The invariance of the Hamiltonian of the color-superconductivity imphes that 

< O >=< VOV^ > (27) 

and 

< O >=< TO^T-^ > . (28) 

The proof of Eqs.|23 ^^'^ H28|l follows from the observation that under or T, a state is either transformed to itself 
( e.g. the ground state ) or another state of the same eigenvalue of H. The hermitian conjugation on RHS of H28fl is 
because of the anti-unitarity of T. See for details. 

The implication of V invariance can be obtained trivially. We have 

S{f,t) ^UpS{-r,t)ul, (29) 

or equivalently 

S-\p,u)^UpS-\-p,Lo)ul, (30) 

which gives rise to 

(l^i-P,^^) = Hp,^)- (31) 

Thus the gap function is even with respect to the spatial momentum. As to the time reversal, Ea. (|28|l together with 
Eqs.lEHll and ^ yield 

Sap{r,t) = I < \T{^x{0,0),¥^ir,t)}T-'\ > {T^)),f3 (32) 
- *(C/t)ap(4)7« < l{*p(0,0),^';(f,-<)}| > (r4)A/3 = mr^)^pSpx{-f,t)iUTri)xf3. 

Therefore 

S{r, t) = C/^r45(-f, i)f7r4 = UlT^S{-r, tprTi, (33) 

which gives rise to the relation 

S-\p,u;) = ulr4S-\~p,Lj)UTT4. (34) 

It follows from 1221) and (EH that 

'/'*(P, -^^) = 0(p,w). (35) 

Consequently, we arrive at 

Theorem 1: The invariance of 2SC under both V and T implies the following off-diagonal structure of the inverse NG 
propagator of real energy ((TH) : 

*(P,^) = -'/*0(p,t^)75A2T2 (36) 
^'(p,^) = ^) 75 A2T2 

where (j){p, oj) an even function of p and satisfies the relation 0(p, uj) — <j)* {p, —lo). 

Decomposing (/)(p, w) into its real and imaginary parts, (f){p,u!) — Re^(p, w) -|- ihn(j)(j), uj) , we find that 
Re(^(p)(Im</)(p, (jj)) is an even(odd) function of tu, a statement that can also be established by analytically contin- 
uating the real solution of the gap equation with Matsubara energy, shown in Ref. [28] and in subsequent sections. 
But the reality of the gap function of Matsubara energy follows from the invariance under V and T as is stated in 
the theorem 2 below. 



B. Matsubara formulation 

The Matsubara quark propagator is defined by 

So.f,{r,T) = [< *„(f,T)*^^(0,0) > eiry < ^'l(0,0)*„(f,r) > di~T)]{Ti)xf, (37) 
The absence of the factor i in comparison with is necessary in order to match the Fourier transformation of l|37|l . 

5,(p)= [\r [ d^re"'^~'P-''S{r,T), (38) 







to the analytic continuation of Eq. ({T^ . i.e. 5y(p) — S{p,w), where i' — 27rT(n + i) is the Matsubara energy. On 
writing 

S^\P)=(^M t^S) (39) 



with ^y{p) = -V*<l>iyi'P}l5hT2 and ^l{p) = il<l>li'P)l5hT2 The self-conjugate relation 

^{f,T) ^ C^{r,t). (40) 

implies that 

S:-\-p)^C^S-\p}C, (41) 

which leads to 

€iP) = €i-P)- (42) 
It follows from (|27|l and the ^'-transformation of the Matsubara operators '^{r, t) and ^{r, t) 

P^'(f, t)V^ = C/p«'(-f, r), (43) 

and 

V^{r, t)V^ = ^{-r, t)uI, (44) 

that 

S-\p) ^ UpS-\-p)uU (45) 

which implies that 

Mp) = M-P)- (46) 
Finally, the T invariance formula H28|l together with the transformation rules 

r«'(r, r)r-i = J7T^'(r, r), (47) 

and 

T^(f, r)r-i = *(f, r)C/^, (48) 

give rise to 

s-\p) = ulrJ-\-p)UTr4. (49) 

It follows then that 

(blip) = <^.(P) (50) 

Combining (|42l (|45|l and H50|) . we end up with 
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Theorem 2: The invariance of 2SC under both V and T imphes the following off-diagonal structure of the inverse NG 
propagator of Matsubara energy : 

= ~V*Mp}15X2T2 (51) 

Kip) = V(t>i'{p)lb^2r2 

where {p) a real and even function of p. 

In the next section, we shall derive the equation for the gap function of real energy, lo) by analytic continuation 
of the gap equation in Matsubara formulation. To unify the notation throughout the continuation process, we shall 
write iv) for (j)^{p). Because of the weak coupling approximation employed, the dependence on p may can ignored, 
leaving the gap a function of energy only. 

The two theorems we have established for 2SC apply to the the gap function of CFL as well. 



III. ELIASHBERG EQUATION FOR GAP PARAMETER AT REAL ENERGY 

The original Eliashberg theory formulated for an electronic superconductor of strong pairing force regards both the 
energy gap and the quasi particle weight ( wave function renormalization ) analytic functions on the complex energy 
plane cut along the real axis. They are determined at equal footing from a pair of self consistent equations of the 
electron self energy with respect to NG basis, one refers to the diagonal part and the other to the off-diagonal part. 
The latter one corresponds to the gap equation in the usual sense. For QCD at asymptotic quark density, however, 
the full complexity of the Eliashberg equations is unnecessary, and one expects the following weak coupling expansion 
of the energy gap function of the real energy w; 



= e-t[/o(.glni^) + g/i (g In + ...] (52) 



for (7 <C 1, but g\n^ = 0{1), where tOc = 0{-^) and powers of Ing are regarded 0(1) in the expansion. Both 
the exponent k and the function fo are known in the literature and are referred to as the leading and subleading 
contributions of the gap function. The leading order imaginary part shows up in the function /i and is therefore 
corresponds to the sub-subleading contribution to the complex gap function. 

To determine the Im/i, we need only one Eliashberg equation, which is the analytic continuation from the gap 
equation of Euclidean energy. For the sake of notation simplicity, we focus mainly on the 2SC case. The generalization 
to CFL case is straightforward and will be addressed at the end of the next section. 

The gap equation of 2SC for either right- or left- handed gap function with Euclidean momentum is given in Ref. 

. , „^ , i — k ■ q 1 + k ■ q (k — q)^ 
+ - Q)( ^ + ^1^)] (53) 

where T/VY^ = TY^ j cfiq/ {2tt)^ in the infinite- volume limit, n labels the Matsubara frequencies — igo = (2n + l)TrT, 

90 n 

and 

Z{qo)^[l + fHMy\qo\')]-' (54) 

is the quark wave-function renormalization factor 1 13lll4lll5j| in the normal phase. It is worth mentioning here that 
this equation does not contain color, flavor and Dirac indices any more, since they have been computed explicitly. 
For the technical details that leads from to see Ref. ^'^s the longitudinal and transverse gluon 

propagators respectively. We will give their expressions later. The feedback of the energy gap to Z{qo), which should 
be determined from the other Eliashberg equation has been neglected. 

To perform the Matsubara sum over quark energies go, we introduce the spectral representations Q. For the gluon 
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propagators 



Ai{p) = -^+ J dTeP°^Ai{T,p} (55) 



l/T 

Atip)^ J dTeP»^AtiT,p) (56) 



A,,t(T,p) = J dijpi^tiiu,p}{[l + nB{uj/T)]e-^^ + neioJ /T)e'^^} (57) 



The spectral densities are given by 

PiA^,p) - pZ'\^,p)6[oj - oJiM] + P?,f (c^,p)e(p - ^) (58) 

where, nsix) = l/(e^ — 1) is the Bose- Einstein distribution function, the expressions of andpf"* will be given 
explicitly later. We also introduce a spectral representation for the quantity^ 

l/T 


oo 

S(t, = y dc^p(w, q){[l - nF{ujlT)]e-'^^ ~ n;^(L^/r)e'^^} (60) 



Where npix) = l/(e^ + 1) is the Fermi-Dirac distribution function. In Ref. 9], Pisarski and Rischke approximated 
the spetral density p(u!,q) by a delta function corresponding to the quasi-particle mass shell, which is appropriate for 
a real gap. In order to extract the imaginary part of the gap, this approximation will not be made here and we shall 
follow the off-shell formulation of the conventional Eliashberg theory that takes the energy as the argument of the 
gapiUl. 

The Matsubara sums over qq can be computed as(p = k — q) 



oo oo 

tJ2MK-QMQ)^ J d£p(e,g){-itanh(^) + J di.pi(i., i3) x [itanh(^)(j-^ 

«o 



1 



1 1 , 1 1 1 1 

+ ,;ootH—){- 



ka — + e ko + v — e 2 2T ko + i^ + e ko — i' — e ko — ly + e 
1 ^ )]} (61) 



OO OO 

T^At(i^-Q)S(g)- j diypt{i^,p) J d£p(£,9)x [ltanh(^)(j-^ 

90 n n 



£ ko — I' — e ko — 1^ + e 



For the first step, we are only interested in zero temperature case. Setting T — leads 

oo oo oo 

tJ2MK-Q)^{Q)= I dvpi{v,p) [ dep{e,q)x{—^-—-- )-l / dep{e,q) (63) 

^ J J ko + V + e ko - V - e p-^ J 
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Tj2MK-Q)EiQ) = J duptiu^p) J de~p{e,q) x {j^^^^ - k^^^^^ ) (64) 

90 

Here we take both the energy and momentum as the arguments of the imaginary part of the gap function. At weak 
couphng hmit, all quark momenta are close to the fermi surface , we set g ~ A; ~ the trivial parts in the Eq. I|53|) 
can be simplified as 

1 + fc-g ^ (fc + - _ ^ 

'i-k-q l + k-q{k-qf _ (P - , . 

2 + 2 (fc_5-)2 4fcg Akqp^ - ^''''^ 

It is now the energy dependence of the spectral densities which provides the interesting phenomena. So the integral 
of q can be written as 

J d^q — J q^ dqd cos 9 dip — J q^dq^^^dpdip ~ J dqpdpdip\k~q=i_i (67) 

The best way to extract the imaginary part of the gap equation, as is shown in Mahan's book, is to integrate q ~ p 
first and then to integrate the polar angles of q. To perform the integral of q — p, we introduce the expression of 
p{e,q). The relation between p{s,q) and is a special case of the Kramers-Kronig relation, when qo approaches 

to the real axis from above. 

p{e,q)^-lm[E{e,q)] (68) 

TT 

For a complex gap function, the spectral density p{s,q) has a finite width. Only when Im0 vanishes, the spectral 
density comes back to the on-shell form that used in Ref.0] 

OO OO OO 

/ dqTy^AiiK - g)S(Q) - - / dupi{v,p) [ d£Ref(e) x (- — ^— ) - ^ / deRef(e) (69) 

J ^ J J ko + !/ + eko-z/-ep^y 



OO OO 

- / dqTy^AtiK - Q)E{Q) - / diyptii^,p) f deRe{{e) x - , ) (70) 

90 n n 



where, /(e) = (/)(e, g)/-\/[e/Z(£)]^ — (\P-{e^q). The minus sign ahead of Ea. l7()|l comes from the contribution of the 
trivial part in Ea. (|66|l . The factor Ref(e) is obtained because of the integral of momenta q. 

j dqlm[E{e, q)] ^ Im[ J dqE{e, q)] = Im[-i7rf(e)] = -7rRef(ir) (71) 

We now perform the analytical continuation, k^ to lu + irj, and take the imaginary part of the Eqs. H69|) . H7U| ) 

Im J dqTY^MK - QMQ) ^ ~signiLj)Tr J deRef{e)pr\H - s,k - q) (72) 



Im J dqTY^MK - Q)E{Q) = sign{u;)7T J deRef{e)pr\H-e,k-q) (73) 
90 n 
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where, we have made use of Eq. (|58(l and ignored the pf"''^ terms, since the two delta functions can not be satisfied 
simultaneously. It means that the pole terms in the spectral densities of gluon propagators give no contribution to the 
imaginary part of the gap function. The expressions of p™* have been given in Ref. j^l in the limit, uj << p « nig. 



pT\u:,k-q) c Ml^fP^ (75) 

■K p^ + [M'^ujy 

where, A/^ = -^m^jm^ = ^^^^ the gluon mass of the HDL propagators for two flavor quark matter at zero 
temperature. Now, we can perform the integral of the polar angles of q 



2^^ 

J pdpd^pr\-,k-<f) - 2M^. J -^^^ ^ I (77) 



The contribution from the longitudinal gluons to the imaginary part can be ignored in comparison with that from 
the transverse gluons if a; << rug, as is required for the validity of the approximations of the gluon propagators. 
Combining Eqs.lO, GSJ ,113) we obtain 



Im0(c^) = sign(c^)^ J deRef(e) (78) 



For the sake of completeness, we also include the real part of the gap function 



where ujc = ^ gs for two quark flavors and ~ ™g- The structure of the gap function is consistent with that in 
0. For instance , the real part is an even function of the energy while its imaginary part is an odd function of the 
energy. 



IV. THE IMAGINARY PART OF THE GAP 

As in the standard Eliashberg theory |^, the imaginary part of the gap function must be derived from two coupled 
equations of Re(/)(tj) and Im(/)(ti;). However, it is important to notice that the forward logaritm in the real part of 
the gap equation does not show up in the imaginary part, which means that Im0 is down by order g relative to Hecf). 
In the leading approximation, we may ignore Imcf) in Eq. H79(l and the RHS of Eq. (|78|l and determine Im0 from 
the approximate solution for Re^ in the literature, i. e. the first term of Eq. (|52|l . For this purpose, we input the 
real part of the gap function for our calculations coming from the results by Schafer and Wilczek Q ■ Although their 
solution is for an imaginary energy , it actually gives the correct, up to the subleading order, results of the real part 
of the gap function(see the Appendix), 

r Ao, if w < Ao 

Re^(w) = <^ (80) 
[ Ao sin[5 In ^] = Aq cos[5 ln(M)], if > Aq 
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FIG. 1: The dependence of lm(j) on lj at /i = 500MeV. The sohd hne shows the approximate solution and the dashed line shows 
the numerical solution. 



with Ao = 2wcexp(-^^ — 9 — 5/(3\/27r). By making use of the Ea. H8()(l . we obtain the analytic leading order 



expressions of the imaginary part of the gap, which is zero for < uj < Aq and is 



gAo-sm[gln— ] 

Z Ao 



(81) 



for Lo > Aq. Here, we have made approximations by ignoring the imaginary part of the gap in the function /(e) in 
Ea. (|78|l , which means the energy qq is on the quasi-particle mass-shell. So the result given by Eq. (|81|l corresponds 
to Eq. (3.177) in^. 

The imaginary part (|78|l follows from a rigorous analytic continuation of the gap equation with Euclidean energy. 
It is instructive to compare this result with the analytic continuation of the approximate gap function in Euclidean 
space, which takes the form 0, Q 



Ao, 



if M < Ao 



Aosin[gln(i^)],if W\ > Ao 



(82) 



The analyticity of the quark propagator for sufficiently small |a;| implies that of (fi^Lu) under the same condition. The 
Euclidean solution (|82|l suggests two logarithmic cuts on the real axis, symmetric with respect to the imaginary axis 
and leaving a gap at the origin. We have then 



(iiy) = Ao sin[§(ln + In 

2 0+ + ii' 



LUc 



0+ - 



-)] 



(83) 

+ iO+ agree exactly with 



for v » Aq. The real part and the imaginary part of its analytic continuation, ii^ - 
Eqs.lHn|l,(IHIl) for UJ » Aq. 

We have also numerically solved the integral Equation (|78() by using the leading order results given by Eq. H80|l. 
The energy dependences of the imaginary part of the gap function at different chemical potentials are depicted in 
Fig. 1 and Fig. 2. We didn't give the results at higher ui because the validity of the approximations for the gluon 
propagator requires that the energy lu can not be too large (see the discussions below the Eq. (|77|l ), and we choose 
the maximum value of lu to be SOOMeV. 

/^From Fig.l and Fig. 2, we know that the contribution from the off-shell behavior of go can be ignored in comparison 
with that from the on-shell behavior, which agrees with the analysis in 2^, i. e. the off-shell behavior only yields 
sufficiently small corrections to Im</). The imaginary part of the on-shell gap function near the Fermi surface was 
calculated in'2^ with a different approach, The prefactor of the result Ea. (|81|l agrees with that in 28], which means 
that the magnitude of the imaginary part is the same as in 28]. 

For comparison. We depict both the approximate analytic result and the numerical result in Fig. 1 and Fig. 2. 
We used the running coupling g{fi) according to the one-loop beta function as used by Schafer and Wilczek 13]. The 
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FIG. 2: The dependence of Im(j> oncu at fi ^ SOOOMeV. 
the numerical solution, respectively. 



The solid line and the dashed line show the approximate solution and 



points of intersection in the abscissa are just the values of Aq. The imaginary part of the gap function is obviously 
a function of lu and is zero for uj smaller than Aq as expected. At extremely large densities, i. e. fi — 5GeV, the 
magnitude of the imaginary part of the gap function is smaller than one MeV. For chemical potential that are of 
physical interest, ^ = 500MeV, the value of the imaginary part can be as large as several MeV. We should caution, 
however, in this regime g is significantly bigger than 1, and higher order corrections are probably important. But to 
the leading order, the numerical results show that the imaginary part of the gap function is down by g relative to the 
real part of the gap function, which is on the order of lOOMeV in weak coupling. 

A positive imaginary part appears odd since one expect naively that the complex quasi-particle pole as given by 
the condition o;^ — (p — /i)^ — A^((jj) = acquires a positive imaginary part for ImA((jj) << ReA(a;), in violation of 
the causality. In this regard, however, one cannot ignore the wave function renormalization Z(uj) even at the leading 
order. The proper condition for the quasi-particle pole to be used reads Z~^a;^ = {p — fJ-)^ + A^(u;), where 



Z{lu) ~ 1 - 5^ 
is the analytic continuation of H54|l We have then 

Imw ~ — TT^Iwl 



In- 



i7rsign(a;) 



Re(j){uj)lm(j){uj) 



with (jj = Recij the solution of uj — sign(w)-\/(p — /i)^ + [Re(a;)]2. On writing x = g\n^ > 0, we find 



Imuj 



2co °\ ^ 



es sin 2x] < 0. 

2 



(84) 



(85) 



(86) 



and the quasi-particle pole is indeed below the physical sheet. It was also suggested in |23| that the imaginary part 
of Z is necessary to shift the quasi-particle pole to the complex energy plane since the pole in his formulation is on 
real axis for Z=real even with nonzero Imt/). We shall come to this difference in the next section. 

Before ending this section, we shall sketch the parallel analysis for a CFL condensate, which is the most favored with 
three quark flavors. Starting with the CFL gap equation with Matsubara energy and following the same procedure 
outlined in section 3, we end up with a pair of gap equation which amounts to replace the function /(e) of 2SC by 



^/i(e) + ^/2(e) 



(87) 



with 



fiie) = 



(j){e,q) 



^[e/Z{£)Y~A4>^{e,q) 
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and 

/.(e) = , ^^^''^ (89) 

The leading order solutions for the real part and the imaginary part remain given by Eqs. H8U|) and (|81l) but with Aq 
replaced by 2^3 Aq. Because of the two gap structure, one may expect three pieces for the solution for three different 
domains, lo < Aq, Aq < uj < 2Ao and 2Ao < uj < 6. But the correction caused by this complication is of higher 
orders. 



V. CONCLUDING REMARKS 



In summary we derive perturbatively the equation of the imaginary part of the gap function for dense QCD in 
weak coupling and solve this equation numerically. The result show that, up to the leading order, the imaginary part 
is small at extremely large densities and can be as large as several MeV for densities of physical interest. We find 
that the imaginary part of the gap function is down by g relative to the real part for the energy < uj << rUg, in 
agreement with the result of Ref. obtained with a different approach. In addition we generalize our result from 
2SC case to CFL case. 

Besides technical differences, the author of (221 made the assumption that the quasi particle energy function e((/i) 
in the function S of l(5^ depends on the imaginary part of the gap according to 

£^(0) = (q - ^)2 + = _ + (Re0)2 ^ (j^^)2^ (90) 

generalizing the treatment of i and where (j) is valued on the quasi-particle mass shell. But he also raised the 
suspicion about this form because it leads to a real quasi-particle pole with a nonzero spetral width 's^. In our 
formulation along the line of Ref. |23| , we have 

€^{4)) = {q- fi)^ + <P^{lu) = (q- ^if + [Re0(a;) + ilriv(t>{u)f (91) 

with (j) valued off shell, since the form |0p cannot be maitained in an analytic continuation. Although the difference 
between eqs. (|9Ul and H91|) does not affect the result of Imcj) up to the order considered in both works, it is conceptual 
and worth clarifying. So we did in Sect. II and we conclude that H91I ) is the right form of the quasi particle energy 
with a complex energy-momentum dependence if the system conserves both V and T. The two NG off-diagonal 
elements are complex conjugate of each other under the same Matsuhara energy only (see Ea. (|42ll ). For real energy, 
the complex conjugation relates the two off-diagonal elements with oppisite sign of energy (see Ea. (|22|l ). In addition, 
the invariance under V and T renders the two off-diagonal elements differ by a constant phase factor only (Theorems 
1 and 2). According to H91|) . the imaginary part of the gap function does shift the quasi-particle pole off real axis. 
But it is still necessary to have ImZ 7^ to place the quasi-particle pole under the branch cut. 

At asymptotically high value of /i, only quarks close to the Fermi surface participate in pairing. The analysis 
carried out here, strictly speaking, applies only to this case with lo << rrig. As fj, is reduced towards the realistic value 
accessible in a compact star, the pairing phase space spread away from the Fermi surface such that the difference 
between uj ^ A and uj ~ nig becomes less clear-cut. The extrapolation of the solution of Ref. [25 and ours to a; ^ 
implies Kecj) ~ Imi/) there. Therefore, the real and imaginary parts of the gap function have to be treated equally and 
one has to stay with Ea. l|91|) for the quasi-particle energy. In addition, the corrections to the spectral density Ht)8|) 
and the correction from the the imaginary part of Z{uj) have to be collected for a more accurate determination of 
lm(j}{u;) [23. 

At realistic value of /i, the pairing channel is speculated to be dominated by instantons and the color superconduc- 
tivity is described by NJL effective action. The magnitude of the gap is expected to be considerably larger than that 
of the one-gluon exchange. While the bare four-quark vertex of NJL gives rise to an entirely real gap parameter, its 
one loop correction contains damping terms and may generate an imaginary part of the gap. In view of the crudeness 
of the weak coupling approximation of NJL, it would be interesting to examine this possibility. 

We expect the nonzero imaginary part of the gap will affect transport properties, the interaction rates, like the 
neutrino emission rates of color superconductors and thus is useful for understanding the evolution and structures of 
a compact star|3lll3ll3ll3ll35l|. 
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APPENDIX A 



In this appendix, we shall solve the gap equation for Re0(aj) up to the subleading order analytically. Neglecting 
the feed back from lm(j){iLj) and the imaginary part of the wave function renormalization|2^., the real part of the gap 
equation is given by 



Re0(c.) = ^ 



de 



Ao 



In- 



In- 



Re<?!)(e) 



subject to the free boundary condition: 



Z-2(Ao)A2 - [Re0(Ao)]2 = 0. 



(Al) 



(A2) 



We shall solve the nonlinear gap equation by iterations starting with a constant gap, Re'"'(/)((jj) = Aq. In each step 
of iterations, we replace Ke4>{uj) inside the square root on RHS of (|A1|) by that obtained from the previous step and 
the integral equation becomes linear subject to the nonlinear boundary condition ljA2|l . To be more specific the linear 
integral equation for the n — th iteration defines an eigenvalue problem 



J An 



(A3) 



for Aq < w, e < Wo of the kernel. 



In- 



(a> + e) 



(A4) 



with the eigenvalue E = 1. The value of Ag is determined upon identifying E with the maximum eigenvalue, i^maxC^o) 
of i. e. 



-Emax(Ao) = 1. 



(A5) 



The normalization of the corresponding eigenfunction, Re0^"''(a;) is fixed by the equation HA2|) and the value of 
Ret/)^"^'^) for < a; < Aq can be obtained from that for w > Aq by means of the integral equation ( IA3I) 

The eigenvalue problem ljA3|l can be analyzed with the perturbation method developed in P, K^. As will be 
justified later, only the first iteration is required for our purpose. The kernel of the gap equation for the first iteration 
reads 



In- 



In- 



Z{e) 



where we have suppressed the superscript oi K{uj,e). Decompose the kernel according to 

K{uj, e) = Ko{uj, e) + K^ico, e) + K^iuj, e) + K^icu, e) + ... 

where 



CJr 1 



with a;> = max(ci;, e) 



ujc Z[e) - 1 



(A6) 



(A7) 



(A8) 



(A9) 
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Ktii^, e) = -g^ln — [ / - -) , (AlO) 

Kli^uj, s) =f\]-( In + In - In ^1 i (All) 

and . . . represents higher order terms. We notice two sources for the logarithm In ^ upon integrating the kernel over 
e with Re(/)(Ao) ^ 0, the one corresponds to the forward singularity of the one-gluon exchange and the one from the 
nonvanishing DOS at the Fermi surface. Both contribute to and render the eigenvalue of the order 0{g^ In^ ■^). 

Then the condition l|A5p implies that In ~ |. Following this rule and Eq. (|54|l . is of the order g. Only the 

forward logarithm contributes to K\ so it is also of the order g. As to K'i, upon a power series expansion of e, both 
logarithms are removed and its contribution is beyond the subleading order. On writing the maximum eigenvalue as 

i^max =Eq+EI+E\^ ... ( A12) 

where Eq is the maximum eigenvalue of Kq and E\ the first order perturbation brought about by ifj^, K\ . We 
have 

i?o = ^ln^^, (A13) 

TT* Ao 

i?? = ^g2ln21n^, (A15) 
According to (|53|l . the gap parameter up to the subleading order reads 

Ao = « . (A16) 

The zeroth order wave function, that solves the eigenvalue problem 

_2 de , LOr, 



eads 



Re(/)(w)=r/ — In— Re0(e), (A17) 
Jao ^ ^> 

Re(/)(tj) = A sin f 5 In =A cos f g In J^V (A18) 
V' \uj\/ \ Aq/ 

with A a constant to be fixed by the nonlinear boundary condition 1A2|I . Since the the first order perturbation of 
the wave function is suppressed by at least by a factor g relative to (|A18|I . we may ignore them for the rest of the 
construction. For the same reason the nonlinear boundary condition can be approximated by Aq = Re0(Ao). It 
follows from the property 

^Re</.(^) (A19) 

that v4 = Ao to the subleading order. Since the RHS of Eq. (|A3(I depends only on Kc(p{uj) for > Ao, we may apply 
it to determine the solution for < w < Aq and the result is Ke(j){uj) ~ Ao. 

The kernel of the gap equation for the next iteration differs from (jA6|) by the dependence of the gap parameter 
inside the square root. Expanding the cosine of (|A18|I . the additional perturbation it brought about is of the order of 

Aq {e^ - Al)2 (e2-Ag)2 
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We have 

r^o i-oo 1 2 

/ deSK{uj,e)^-9' dx (A21) 

JAo J I (x'' - 1)2 

Therefore its contribution to the eigenvalue is beyond the subleading order. 
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